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Математическогєе моделированиеє оптимизации системьі 
«деготь - полимер - древесньй гидролизньтай ЛИГНИН» 


У статті подана розроблена математична модель оптимізації системи «дьоготь -- полімер (полівінілхлорид, полі- 
стирол) -- деревний гідролізний лігнін». Здобуто оптимальні концентраційні співвідношення між компонентами, 
що складають систему, за допомогою методів математичного планування експерименту. 
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В статье изложена разработанная математическая модель оптимизации системь «деготь -- полимер (поли- 
винилхлорид, полистирол) -- древесньгй гидролизньй лигнин». Полученьт оптимальньюе концентрационньте со- 
отношения между составляющими систему компонентами с помощью методов математического планирова- 
ния зксперимента. 

Ключевьге слова: модель математическая, параметрьт оптимизации, функция отклика. 


Кам'яновугільні дьогті є слабко концентрованими суспензіями (а- 1 В-фракція в 
суміші низькомолекулярних ароматичних вуглеводнів), для яких характерна невисока 
в'язкість, низька когезія, відсутність еластичності, схильність до інтенсивного старіння | 1. 
Ефективним заходом поліпшення якості є їх фізико-хімічна модифікація полімерами, а 
також одночасна дія на дисперсну фазу і дисперсійне середовище органічних в'яжучих 
коксохімічного виробництва комплексними добавками, що містять полімери (полі- 
вінілхлорид, полістирол), які суміщаються з ними, та активні дисперсні наповнювачі 


21, 131. 
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Як параметри оптимізації дьогтеполівінілхлоридної системи |4) було прийнято: 

- температуру розм'якшення (не менше 34"С); 

- температуру крихкості (не вище 13"С); 

- інтервал пластичності (не менше 47"С); 

- еластичність (понад 45"С). 

Регресійний аналіз виконано за програмою мовою МаїСАР з використанням 
пакета 5 ВЕЕК Юг У/ппдоууз. Аналіз діаграм дозволив апроксимувати криву температу- 
ри розм'якшення поліномом другого ступеня, а криві температури крихкості 1 
еластичності описати неповним поліномом третього ступеня (У, ХУ», ХУ з). Здобуто 
такі Р регресії (41: 

- 32,51. 6,115х1 1 7,796-х» - 0,266-ху:х» Ч 0, 1609-хі; - 1,331: ху 

й чо -15,34 Я 3,094: ХІ а 1,978: Х2 ка 0,912: ХХ к 1,156: Ж Я 1,57: УЮ Я 1,031:хухо 

У з з 68,19 - 17,0хху:х» - 19,34-сху - 11,47:ха" - 16,98-ху "ха - 18,19:х" - 8,67:хо". 

Встановлено, що оптимальною консистенцією середовища, яке модифікується, є в'яз- 
кість СУ - 150-200 с, аоптимальна концентрація полівінілхлориду становить 1,5 -- 2,0 9/6 |41. 


Суміщати полістирол з кам'яновугільним дьогтем слід при температурі 105 -- 110"С 
впродовж 40 - 50 хвилин |5|. Оптимальною в'язкістю вихідної сировини є в'язкість 
СУ -10 - 30 с, араціональний вміст полістиролу має бути 0,5 -- 2,0 9) 6). 

Рідкі кам'яновугільні дьогті, у яких енергія теплового руху дисперсійного сере- 
довища (у - фракція) дорівнює або більше енергії зв'язку у вузлах макромолеку- 
лярної сітки, що утворена за рахунок механічного зачеплення 1 переплетіння над- 
молекулярних структур полістиролу (неполярного полімеру), слід зміцнювати ліофіль- 
ними наповнювачами, наприклад, пеком, сажею, молотими довгополуменевим вугіллям, 
деревним гідролізним лігніном, золою-виносу ТЕЦ або окислювати |7). 

В роботі |8)| доведено, що деревний гідролізний лігнін має високу здатність, що 
структурує, дьогтеполімерних в'яжучих. На комплексні кам'яновугільні композиції впли- 
вають такі фактори: в'язкість вихідного дьогтю; концентрація полімеру; концентрація 
наповнювача. 

Комплексні в'яжучі є складними сполуками. Тому виникає необхідність їх 
оптимізації за допомогою методів математичного планування експерименту зі здобуттям 
необхідних їх фізико - механічних властивостей. 

В даній роботі як об'єкт дослідження було прийнято: 

- середовище, що модифікується, - кам'яновугільні дьогті - складені з 
середньо температурного пеку і антраценового масла; 

-- полімери - первинні відходи виробництва полівінілхлориду (ПВХ) Дніпродзержинсь- 
кого об'єднання «Азот» з молекулярною масою 12-10" в.о і розміром часток (6,3 -- 63)10"м і 
полістиролу (ПО) ВАТ «Концерн «Стирол» (м. Горлівка Донецької області) з молекуляр- 
ною масою 9-10" в.о і частками розміром менше 6,3-10"м. 

- наповнювач - деревний гідролізний лігнін (ДГЛ) Бобруйського гідролізного 
заводу (Білорусь) - відхід спиртового 1 дріжджового виробництва, який одержують у 
вигляді осаду хвойних и листяних порід деревини методом гідролізу (розбавленою 
сірчаною кислотою). 

Комплексне кам'яновугільне в'яжуче готували суміщенням кам'яновугільних дьогтів з 
полімером при температурі 115 - 125"С впродовж 30 хвилин |8|. Потім додавали 
порошкоподібний деревний гідролізний лігнін і ще перемішували 30 хвилин. 

Для оптимізації системи «дьоготь - полімер - деревний гідролізний лігнін» 
було використано композиційний несиметричний план |9| на трьох цілочисельних 
рівнях (- 1; 0; --1) з коефіцієнтом кореляції гі, 2 0,10, де 1,) - 1,23. 

Оптимальні склади систем «дьоготь - ПВХ - ДГЛ» 1 «дьоготь - ПС - ДГЛ» ви- 
значали як оптимальні області допустимих значень факторів Х., Х», Ха (табл. 1). 
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Таблиця 1 - Значення варійованих факторів 
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Фізичний зміст фактора 
Мо Умовна | Масова Масова 
з.п. Система в'язкість | концентрація | концентрація 
дьогтю полімеру, Ус | наповнювача, 
по С,с о 
Х Х, Х; 
1 | Дьоготь-ПВХ- Інтервал варіювання 100 1 20 
ДГЛ Рівні фактора -1 50 0 0 
0 150 1 20 
ЧІ 250 2 40 
2 | Дьоготь-ПС-ДГЛ | Інтервал варіювання 100 245 20 
Рівні фактора -1 50 0 0 
0 150 259 20 
ЧІ 250 5 40 


Ці оптимальні області обмежені поверхнями рівня функції відклику з 
кожного з параметрів оптимізації (табл. 2). 


Таблиця 2 - Параметри оптимізації системи 


Хо зп. Код Фізичний зміст параметра оптимізації | Одиниці | Граничні зна- 
параметра виміру | чення функції 
оптимізації відклику 

1 У! Оптимальний вміст в'яжучого в суміші У Не більше 8,5 
(в перерахунку на дьоготь) 
2 У» Температура розм'якшення "Є Не менше 33 
3 У Еластичність при 0"С о Не менше 30 
4 У Границя міцності при стиску при 20"С МПа  |Неменше2,5 
5 У Границя міцності при стиску при 50"С МПа  |Неменше 1,0 
6 У є Границя міцності при стиску при 0"С МПа  |Небільше 12,0 
7 У з Коефіцієнт водостійкості - Не менше 0,8 
при довготривалому водонасиченні 
8 Уз Відношення границі зсувової міцності - 11-11, 
до напруження сталої течії 
Матриця планування експерименту подана в табл. 3. 
Таблиця 3 - Матриця планування 
Мо з.п. Х Х, Х. 
1 -1 -1 -1 
2 0 -1 0 
З 1 -1 1 
4 -1 0 0 
5 0 0 1 
6 1 0 -1 
7 -1 1 1 
8 0 1 -1 
9 1 1 0 
10 -1 -1 1 
1 1 -1 -1 
12 -1 1 -1 
13 1 1 
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Загальний вид функції відклику являє собою рівняння другого порядку. Це рівняння 
характеризується такими статистичними параметрами: дисперсія адекватності 5 - 0,016; 
коефіцієнт варіації 6 - 2,3490; кореляційне відношення КВ - 0,984. 

Здобуто такі рівняння регресії: 


Система «дьоготь - ПВХ - ДГЛ» 


У - 8,239 -- 0,217:Х, -- 0,094- Х» -- 1,992: Х3-- 0,076: Ху: Х» -- 0,055: Ху: Хз- 
0,071: Ху Ха 0,133.Ху чн 0,020-Х»7 -- 0,986: Хо 

У» - 28,335 4 3,812: Х; -- 4,931: Х» -- 6,950. Ха - 1,501: Ху: Х» - 0,309: Ху: Хз- 
2,115. Ху Ха - 0,674-Х у" - 0,079-Х,7 - 1,453: Ха 

УЗ з 13,572 4 3,016-Х,; - 24,575- Х» - 10,111: Ха -- 1,112. Ху: Х» - 1,755: Ху: Хз- 
8,160 Хо Ха 4,394-Х 2 -- 1,192:Х52 4 7,870. Хо 

Уд - 2,751. 0,311:Х; я 0,693: Х» - 0,592. Ха - 0,122: Ху: Х» З 0,154: Ху: Хо 
0,186 Х»- Ха - 0,113:Х7 4 0,098.Х»7 - 0,131: Хз 

У 5 -- 0,938 -- 0,060-Х, -- 0,274. Х» -- 0,158: Ха -н 0,001: Ху Х» - 0,016: Ху Ха 
0,102: Х»- Ха - 0,012:Х,7 - 0,343:Х» - 0,050: Хз; 
У в -- 9,061 -- 1,367: Хі - 1,838: Х» -- 2,324: Ха - 0,343: Ху: Х» чн 0,524: Ху Ха 
0,326: Х»: Ха - 0,611:Х7 - 0,460:Х». - 0,759. Х5.; 
У з - 0,809 -- 0,032:Х; -- 0,085. Х» -- 0,023: Хз - 0,024: Ху: Х» -н 0,006: Ху: Ха 
0,001 Х»- Ха - 0,027:Х,7 - 0,349-Х»- - 0,023: Хз 

У в з 11,258 -- 0,0:Х; -- 0,108: Х» -н 0,157: Ха - 0,033 Ху: Х2 3 0,04 Ху Хо 
0,094 Х» Хз - 0,133Хи. - 00-Х» - 0,0. Хз.. 


Система «дьоготь - ПС - ДГЛ» 


У - 8,062 -- 0,161:Х, - 0,051 Х» -- 2,087: Ха - 0,031: Ху: Х» - 0,012: Ху: Хз- 
0,007 Х»- Хан 0,176. Ху. чн 0,216-Х»7 0,738: Хо 

У» «з 33,801 -- 4,413:Х, -- 4,497. Х» 4 7,634: Ха - 1,343: Ху: Х» - 0,966: Ху: Хз- 
1,133 Ху Ха - 4,973:Ху/ -н 0,443:Х» - 5,211: Ха 

У з з: 33,831 4. 0,462:Х, 4 21,069- Х» - 1,093: Ха -- 0,929- Ху: Х» - 1,691: Ху: Хз- 
0,620: Х»" Ха 1,876-Хиу. - 11,088:Х»- - 3,969- Ху; 

уд з 2,736 4 0,383: Х - 0,570. Х» - 0,486: Ха 0,056: Ху: Х» -- 0,096- Ху: Хо в 
0,066 Х» Хз - 0,255. Ху. -- 0,079-Х»/ - 0,101: Ху; 

у 5 -- 0,918 -- 0,092-Х, -- 0,198. Х» -- 0,149. Ха -н 0,033: Ху: Х» - 0,059: Ху: Ха 
0,076 Х» Хз - 0,016. Ху. - 0,124-Х»7 - 0,012: Хз. 

Ув - 9,325 -- 1,538: Хи 2,139. Х» -- 1,723: Ха - 0,277: Ху: Хо - 0,108: Ху: Хз- 
0,180 Х»- Ха - 1,365: Х. - 0,459-Х» -- 0,031: Хз 

У з - 0,826 -- 0,016-Х; -- 0,086: Х» -- 0,051: Хз - 0,040: Ху Х» - 0,006: Ху Хо з 
0,025. Х»" Ха - 0,08: Хи. - 0,045-Х»/ - 0,002: Ху. 


Рівняння регресії перевірено на адекватність, 1 вони відповідають критерію Фішера. 

Для визначення спрямування обмеження для кожної функції відклику в допусти- 
мій області було прийнято значення, близьке до граничного. За програмою для 
зазначених двох величин кожного з параметрів оптимізації У; - Уз було обчислено 
значення факторів Х, Х», Хз, які задовольняють рівнянням функції відклику для 
кожної системи. Послідовно фіксували один з факторів на рівнях: -Ї; -0,5; 0; 1,0. 
Іншому фактору було надано фіксовані значення з тим самим інтервалом і було 
визначено третього фактора як корені рівнянь, що задовольняють функції відклику. 
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Рисунок 1 -- Площини рівнів систем: 
а)дьоготь С19 - 50с- ПВХ - ДГЛ; 
б)дьоготь С35 - 100с- ПВХ ДГЛ; 
в)дьоготь С15 - 150с- ПВХ ДГЛ; 
г)дьоготь С14 -- 200с- ПВХ - ДГЛ; 
д)дьоготь С10 - 250с- ПВХ - ДГЛ 
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Рисунок 2.-- Площини рівнів систем: 
а) дьоготь С15 - 50с- ПС-- ДГЛ; 
б) дьоготь С19 - 100с- ПС - ДГЛ; 
в) дьоготь 014 - 150с- ПС - ДГЛ; 
г) дьоготь С35 - 200с- ПС - ДГЛ; 
д) дьоготь С15 з 250с- ПС - ДГЛ 
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Масова концентрація , ПВХ, 90 


в) 
Рисунок 3 - Діаграма оптимальний вміст в'яжучого в суміші 
(в перерахунку на дьоготь) -: 
а) умовна в'язкість дьогтю по С. і масова концентрація ПВХ; 


б) умовна в'язкість дьогтю по С. і масова концентрація ДГЛ; 
в) масові концентрації ПВХ і ДГЛ. 
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6) 


Температура 
розм'якшен- 
ня, "С 


Масова 
-- концентрація 


дл, 


Масова концентрація , ПВХ, 90 
в) 
Рисунок 4 - Діаграма температура розм'якшення - : 
10 


а) умовна в'язкість дьогтю по С і масова концентрація ПВХ; 


б) умовна в'язкість дьогтю по С, і масова концентрація ДГЛ; 


в) масові концентрації ПВХ 1 ДГЛ. 

Рівняння У - (Х., Х», Хз) функції відклику, що пов'язує варійовані Х., Х2, Хз, за 
суттю є рівняння гіперболічної поверхні другого порядку в чотиримірному просторі У, 
Х, Х», Хз. Фіксуючи значення У, одержуємо поверхню рівня функції відклику в просторі 
факторів Х., Х», Хз. Потім, фіксуючи значення фактора, наприклад, Х/, маємо лінії рівня 
функції відклику на площині варійованих Х», Хз. 
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Порівнюючи між собою положення ліній рівня для граничного і приграничного 
допустимих значень параметра оптимізації, було визначено допустимі області значень 
факторів за певним параметром оптимізації. Потім було зведено в одну площину рівня 
орієнтовані лінії рівня кожного з параметрів оптимізації, що дає можливість виділити в 

Для параметрів оптимізації, що обмежують допустиму область, всі п'ять площин 
рівня з побудованими лініями рівня було об'єднано в плановій аксонометричній проекції. 
Порівнюючи між собою одержані поверхні рівня параметрів оптимізації, що обмежують 
допустиму область, було побудовано лінії перетину цих поверхонь між собою, а потім у 
вигляді об'єму було виділено область допустимих значень факторів, які оптимізують 
склад системи. Встановлено, що для системи «дьоготь - ПС - ДГЛ» обмежуючими 
параметрами оптимізації є У, У», У з та У (рис. 2), а для системи «дьоготь - ПВХ - ДГЛ» - У! 
та У» (рис. 1). Аналіз діаграм функцій відклику, наприклад, У! (рис. 3) та У. (рис. 4) 
доводить, що оптимальний вміст в'яжучого в суміші майже не залежить від умовної 
в'язкості дьогтю (рис. 3, а), але значно залежить від концентрації полівінілхлориду 
(рис. 3,а, в). В разі низької умовної в'язкості дьогтю деревний гідролізний лігнін 
суттєво впливає на оптимальний вміст в'яжучого в суміші і майже не вливає на 
нього для підвищеної в'язкості дьогтю (рис. 3, б). 

Зростання умовної в'язкості дьогтю і концентрації полівінілхлориду призводить до 
збільшення температури розм'якшення рівнозначно (рис. 4, а). Підвищення концентрації 
деревного гідролізного лігніну в малов'язких дьогтях майже не впливає на температуру 
розм'якшення, а для дьогтю підвищеної в'язкості спостерігається її помітне підвищення 
(рис. 4, б, в). В системі «дьоготь - ПВХ - ДГЛ» лише при умовній в'язкості С - 130-250 с 
(рис. 1) спостерігаються оптимальні області концентрації полімеру і наповнювача при 
граничних значеннях прийнятих факторів. Вміст полімеру в дьогті понад 290 за масою 
призводить до перенасичення розчину (дьоготь -- ПВХ). Набряклі частки ПВХ, які ви- 
діляються з розчину, є додатковими вузлами зшивання просторової сітки, що обумовлює 
зміцнення системи. При цьому, чим менше в'язкість модифікованого дьогтю, тим 
інтенсивніший тепловий рух і, як наслідок, для зміцнення структури в'яжучого необхідна 
більша кількість полімеру. Масова концентрація деревного гідролізного лігніну для всіх 
систем приблизно однакова (2096) і обмежується таким параметром оптимізації, як 
оптимальний вміст в'яжучого в суміші (У |). 

Оптимальні ж області системи «дьоготь -- ПС - ДГЛ» дещо розвинуті і виявляють- 
ся при в'язкості дьогтю С. - 150 та С - 200 с. В цих системах контакти між частками 


наповнювача здійснюються за допомогою надмолекулярних структур неполярного 
полімеру (полістиролу), концентрація якого в рідких системах досягає 5906 за масою. 
Таким чином, оптимальними складами досліджуваних систем є: 
Г) дьоготь в'язкістю С. - 130 - 250 с - 10090; ПВХ - 1,6-- 2,090; ДГЛ - 13-219б; 


2) дьоготь в'язкістю С! - 125 - 225 с - 10090; ПС - 2,5 - 4,490; ДГЛ - 16-2290. 
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